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Abstract. Electrical and optical properties of binary inhomogeneous media are currently 
modelled by a random network of metallic bonds (conductance uo, concentration p) and 
dielectric bonds (conductance cti , concentration 1 — p). The macroscopic conductivity of 
this model is analytic in the complex plane of the dimensionless ratio h = ai/ao of the 
conductances of both phases, cut along the negative real axis. This cut originates in the 
accumulation of the resonances of clusters with any size and shape. We demonstrate that 
the dielectric response of an isolated cluster, or a finite set of clusters, is characterised 
by a finite spectrum of resonances, occurring at well-defined negative real values of h, 
and we define the cross-section which gives a measure of the strength of each resonance. 
These resonances show up as narrow peaks with Lorentzian line shapes, e.g. in the weak- 
dissipation regime of the RL — C model. The resonance frequencies and the corresponding 
cross-sections only depend on the underlying lattice, on the geometry of the clusters, 
and on their relative positions. Our approach allows an exact determination of these 
characteristics. It is applied to several examples of clusters drawn on the square lattice. 
Scaling laws are derived analytically, and checked numerically, for the resonance spectra 
of linear clusters, of lattice animals, and of several examples of self-similar fractals. 
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1. INTRODUCTION 

Frequency-dependent electrical and optical properties of inhomogeneous media are 
often modelled by networks of random complex impedances. The case of binary composite 
media has attracted a lot of attention (see ref. [1] for a review). Each bond of a regular 
lattice is independently attributed a random complex, frequency-dependent conductance 
(or admittance) according to the binary law 

{(Jo with prob. p, 
(Ti with prob. 1 — p 

The dimensionless complex ratio of the conductances of both phases, 

h = 

and the concentration p are the essential parameters of the model. 

The ^ 1 regime describes several situations of interest. As far as static (DC) prop- 
erties are concerned, the h — limit embraces the two well-known cases of the conductor- 
insulator mixture (ui = 0), and the superconductor-conductor mixture (uo = oo). In both 
limiting situations the macroscopic conductivity E exhibits critical behaviour around the 
percolation threshold pc, 

S(c7i = 0) ~ CTo(p-Pc)* {P^Pt)^ 
S((7o = oo) ~ (7i(pc -P)~* (P^Pj), 

where s and t are universal critical exponents. 

Frequency- dependent (AC) electrical and optical properties of metal-dielectric mix- 
tures have attracted much interest recently. Two models have been mostly investigated, 
the R — C model and the RL — C model [1] . The dielectric component is modelled by 
perfect capacitors, while the metallic bonds consist of a resistance in the first case, of an 
inductance in series with a resistance in the second case. The RL — C model has been 
proposed to describe optical properties of several kinds of inhomogeneous materials, such 
as cermets [2] . These materials undergo an optical transition at some value p* > Pc of the 
metallic fraction, where their dielectric constant changes over from inductive to capacitive. 
In both R — C and RL — C models, the frequency dependence of the ratio h is such that 
the low-frequency regime corresponds to \h\ <^ 1 (see section 2.2). 

In the critical region, namely for small h and near the percolation threshold, the 
conductivity obeys a universal scaling law of the form 

^^ao\p-pc\''^±[h\p-pc\-'^'+'^) (|/i|«l, |p-Pc|«l), (1-4) 

where ± refers to the sign of {p — Pc)- 

The present work aims at a better knowledge of the analytic structure of the con- 
ductivity in the complex /i-plane, for any fixed concentration p. In spite of its apparent 
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(1.1) 



(1.2) 



(1.3) 



academic character, this point is essential for a quantitative understanding of both the 

transient response of the R — C model to a time-dependent excitation, and the resonance 
spectrum of the RL — C model [1]. Hereafter we shall mostly focus our attention on the 
second situation. 

Consider first a finite network, defined by putting on the bonds of a graph conduc- 
tances which assume two values, do or ai. The conductance of this network between any 
two nodes assumes the form Y = aoF{h), where F{h) = N{h)/D{h) is a rational function 
of the complex ratio /i, i.e., N{h) and D{h) are polynomials, whose degrees are roughly 
equal to the number of bonds of the network. Furthermore the zeros of the conductance 
(i.e., the zeros of N) and its poles (i.e., the zeros of D) are negative real numbers, and 
they alternate, namely there is exactly one zero between any two consecutive poles, and 
vice-versa. These properties are at the origin of various analytic representations and rig- 
orous inequalities for the conductivity and the dielectric constant of random media. This 
body of knowledge is currently referred to as the Bergman-Milton theory [3] . 

As a consequence, the singularities of the macroscopic conductivity take place for h 
real negative, so that E is analytic in the complex /i-plane cut along the negative real axis. 
We introduce the notation 

Disc E(/i) = I (p{h + zO) - E(/i - zO)) = 2 Im E(/i + zO) (1.5) 

for the discontinuity of the conductivity along this cut. 

It is interesting to first look at the prediction of the effective-medium approximation 
(EM A). This old and very commonly used approximate scheme [4-6] amounts to resum- 
ming the one-impurity effects in a self-consistent way. In the case of the square lattice the 
EMA formula reads 

S(EMA) ^ ^^(^(^ _ 1/2) (1 - /i) + 1/2)2(1 -/i)2 + /i). (1.6) 

This expression shows that the conductivity is cut along a finite interval [/imin, /imax], with 

1 + 4p(l -p)±4. ^p{l-p) 
{2p - 1) 

and that its discontinuity reads 

DiscE^E^^) = 2ao\p - 1/2| V(/imax - /i)(/i - /imin) (/imin <h< h^^^). (1.8) 

It has been argued [1] that the exact conductivity of the random binary model has a 
cut with a non-vanishing discontinuity along the whole negative real axis. The bulk of the 
discontinuity is expected to be roughly given by eq. (1.8), while its tails are expected to be 
much weaker, and to become exponentially small as either /i^O^ or h ^ —oo, in analogy 
with the Lifshitz tails of the density of states of electrons and phonons in disordered solids 
[7]. The following law for D-dimensional lattice problems has been proposed in ref. [1]: 

DiscE - exp (^-Cl/if-^/^) (h^Q-), (1.9) 
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while a more recent investigation [8] rather suggests the above law with formally D = 1, 
independently of the dimension D. 

We propose to get a novel kind of physical insight into the cut of the conductivity, by 
viewing it as a result of the accumulation of the resonances of clusters or sets of clusters 
of any size and shape. The analysis of the complex conductivity of the RL — C model 
in terms of resonances has been tackled in refs. [1, 8-10]. The present work provides a 
quantitative analysis of the dielectric resonance spectrum of any cluster drawn on a regular 
lattice. The case of the square lattice is considered for definiteness. The setup of this paper 
is as follows. Section 2 contains general results on the resonance frequencies of arbitrary 
finite clusters and sets of clusters, drawn on the square lattice, as well as the associated 
resonance cross-sections. Several applications are then presented in section 3, including 
two coupled bonds, linear clusters, arbitrary clusters of a given size (lattice animals), and 
examples of self-similar fractal clusters. Section 4 contains a short discussion. 

2. GENERAL RESULTS 
2.1. Definitions 

The aim of this section is to present a general approach to the resonant dielectric 
response of clusters embedded in a regular lattice. In the following a cluster is any finite 
connected set of bonds drawn on a lattice, and a set of clusters is a finite collection of 
such clusters. We denote by tib the total number of bonds (links) contained in a cluster 
or in a set of clusters, and ns its total number of sites (vertices, or nodes) . For the sake of 
simplicity, we restrict ourselves to the two-dimensional square lattice, spanned by the unit 
vectors {ei, 62}. We take the lattice spacing, i.e., the length of the bonds, as our length 
unit. 

In order to study the dielectric response of a set of clusters, we consider the geometry, 
shown in Figure 1, of a rectangular sample of size M x A^, between two straight parallel 
electrodes. The bonds which belong to the set of clusters, shown as thick lines and referred 
to as impurity bonds, are assumed to have a complex, frequency-dependent conductance (or 
admittance) cri, while the other bonds of the lattice, shown as thin lines and referred to as 
matrix bonds, have a different conductance (Tq. The complex ratio h of both conductances, 
defined in eq. (1.2), will play a central role in the following. 

Our starting point is the Kirchhoff equation for the electric potentials. 



We employ the following notations: x = xiBi + X2G2 = (^ci, X2) is any site of the lattice, 
and y(x) denote the four nearest neighbouring sites of x, i.e., y = x ± e^, with |U = 1 or 
2. The electric potential at site x is denoted by Vy^, and ax,y = CTy^x is the conductance of 
the bond joining the neighbouring sites x and y. Finally the source term Xx is the current 
flowing from the generator into the network through node x; it is non- vanishing only when 
X belongs to either electrode. 




(2.1) 



y(x) 
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Both and Xx have to be solved from eq. (2.1), with the boundary conditions Vx, = 

on the left electrode (x G £"0); and V^^ = V = MS on the right one (x G with S being 
the uniform electric field applied to the sample. The complex conductance (or admittance) 
between the electrodes reads 

where I is the total current across the sample, 

X = - 5^ Xx = 5^ Xx. (2.3) 



2.2. The RL-C model 

One situation of special interest is the RL — C model of inductive (and weakly dissi- 
pative) clusters in a dielectric matrix, already mentioned in section 1. The bonds of the 
clusters consist of an inductance L in series with a weak resistance -R, while those of the 
matrix are perfect capacitances C. The complex conductances at frequency / = uj/{2tt) 
thus read ^ 

o-Q = iCu, ai = ry , T ■ (2-4) 

R + iLu) 

Along the lines of refs. [1, 2, 9, 10], we introduce the microscopic resonance frequency 
(plasmon frequency) 

0.0 = (2.5) 



and the quality factor 



which is a dimensionless measure of the dissipation rate. We also introduce the reduced 
frequency 

y=-, (2.7) 

(jjQ 

so that ^ ^ 

^ " -LCa;2 + ^^cuj ^ -y^ + iy/Q' ^^'^^ 

In the following, we shall mostly consider the regime of a weak dissipation, corre- 
sponding to a large quality factor. In this regime, we have therefore 

^^-^-TTs (<3»1)- (2.9) 

We have recalled in section 1 that Y{h) is a rational function, with alternating poles 
and zeros along the negative real axis. Since the variable h as given by eq. (2.9) is very 
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close to this negative real axis for Q ^ 1, we can anticipate that the poles and zeros lying 
there will strongly affect the frequency-dependent response of RL — C clusters, in the form 
of narrow resonances. This is what we shall now demonstrate explicitly. 

2.3. Resonance frequencies 

We define the resonances of a cluster, or a set of clusters, as the values of h such that 
the conductance Y{h) vanishes, in the limit of an infinitely large network. This definition 
will be justified in section 2.4. 

At a resonance, eq. (2.1) must have a non-trivial solution localised around the 
clusters, in the absence of sources. This equation with 2x = can be recast as 

-(AF)x = (l-M (^x-VV). (2.10) 

yec(x) 

The notation y e C(x) means that the bond (x, y) belongs to the set of clusters, and A 
denotes the finite-difference Laplace operator on the square lattice, defined as 

-{AV)^^Y.{y^-Vy). (2.11) 
y(x) 



The solution of eq. (2.10) with appropriate decay properties at infinity reads 

^^-=E E G^AVy-V^)^ (2-12) 

yeCzeC(y) 

where we have set 

A=^, (2.13) 

and where Gx,y = G(x — y) is the Green's function of the Laplace operator on the infinite 
square lattice. Its main properties are recalled in the Appendix. 

By expressing the consistency of eq. (2.12) for x being a site of the cluster set, we 
arrive to the following characterisation of the resonances: A has to be an eigenvalue of the 
square matrix M, of size ns x ns, defined by 

Mx,y= E (Gx,y-Gx,z). (2.14) 
zGC(y) 

This matrix is not symmetric. It can nevertheless be recast, using a bond representation, 
in the form of a real symmetric matrix of size ub x ns, whose spectrum is manifestly real 
[10]. 

The ns eigenvalues of the matrix M lie in the range < A < 1. Only those different 
from the endpoints (A 7^ and 1) correspond to physical resonances. Wc denote by riR 
the number of these non-trivial eigenvalues, which we assume to be ordered as 

< Ai < A2 < • • • < A„^ < 1. (2.15) 



6 



Even for a single cluster, there is no simple relationship between the numbers of sites, of 
bonds, and of resonances, apart from the inequalities ri]^ tll < ns- Furthermore, some 
of the eigenvalues may be degenerate, i.e., occur with a non-trivial multiplicity. 

The spectrum of dielectric resonances A^ (a = 1, ■ ■ ■ , n^) of a given set of clusters 
thus only depends on the shape of the clusters, and of their relative positions and orienta- 
tions. In the weak-dissipation regime of the RL — C model, introduced in section 2.2, the 
resonances occur at well-defined resonance frequencies fa = a;o/(27r), given by 

Aa = :r^ = ^|y = i.e., coa^u;oJy\- (Q»l). (2.16) 

2.4. Resonance cross-sections 

We now turn to the determination of the analytic form of the conductance Y{h), in 
the regime of a weak dissipation {Q 3> 1). We restrict the analysis to the case where 
the linear sizes M and N of the sample are much larger than the diameter of the set of 
clusters under consideration, so that the resonances are very close to those determined in 
section 2.3, corresponding to clusters embedded in an infinite lattice. Section 4 contains a 
qualitative discussion of finite-size corrections. 

Let M be the matrix associated with the clusters, and let Aa (a = 1, ■ ■ ■ , ur) be its 
non-trivial eigenvalues. We assume, for the sake of simplicity, that all eigenvalues are 
non-degenerate. For each eigenvalue A^, we denote by ha, i/a, and coa the corresponding 
values of the various variables defined above, the last two pertaining to the RL — C model 
[see eq. (2.16)]. We also introduce the associated left and right eigenvectors La and Ra, 
with components x and -Ra,x in the site representation (2.14). These eigenvectors are 
supposed to be normalised in such a way that 

La-Rb = ^ La,^Ra,^^ Sa,b (a, 6 = 1, . • ■ , n/j), (2.17) 

with 6a,b being the Kronecker symbol. 

We look for a solution of the Kirchhoff equation (2.1) in the form 

^x = £xi + W^, (2.18) 

where the first term is the potential in the absence of the set of clusters, with xi being 
the co-ordinate of the node x along the applied field S, and is a perturbation of the 
potential, localised around the clusters, and vanishing on the electrodes. The Kirchhoff 
equations for the potential Wj^ on the clusters can be recast as 

AWx - Yl Mx.yW^y =^Y1 Mx,y2/1. (2.19) 
y€C yec 
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Whenever A comes close to one of the eigenvalues A^ of the matrix M, the solution of eq. 
(2.19) diverges as 

^ A{X)Ra,^, (2.20) 

with Ra being the corresponding right eigenvector. The divergent prefactor A{X) is then 
determined by multiplying the corresponding left eigenvector La through eq. (2.19). We 
thus obtain 

A{X)^-^SCa (A^A„), (2.21) 

with 

jCa=Yl (2.22) 

xec 

This quantity is independent of the position of the set of clusters inside the sample, owing 
to the identity ^x€C -^a,x = 0. The total current across the sample is then evaluated as 

I=ao\N£+^^^ J2 {xi-yi){W^-Wy)\. (2.23) 

(x,y)ec 



The conductance of the sample in the regime of a weak dissipation can now be esti- 
mated by superposing the independent contributions of all resonances to the intensity. We 
are thus left with the following formula 



with 



and 



7a = 



Xail-XaV 



T^a= 5Z (^1 ~ 2/l)(-^a,x - -Ra.y)- 

(x,y)ec 



(2.25) 
(2.26) 



Eq. (2.24) is our main result. The prefactor {N/M)ao is the conductance of the 
matrix sample, in the absence of the clusters. To this background response are superposed 
the ur resonances of the set of clusters. In agreement with the general properties recalled 
in section 1, each resonance shows up as a doublet, consisting of a pole situated at h = ha 
and a zero situated at /i ~ /ia(l — 7^/ (MA'")). 

The strength of each resonance is measured by the distance between the pole and 
the zero, i.e., by the residue 'jaha/{MN). It is inversely proportional to the area MN of 
the sample, and proportional to the parameter 7„, that we call the cross-section of the 
resonance. The latter is indeed interpreted as the area of the region around the clusters 
where the perturbation W^^ of the potential at resonance takes appreciable values. 
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The general result (2.24) can be made more explicit in the Q ^ 1 regime of the 
RL — C model, introduced in section 2.2. In this situation, it is advantageous to consider 
the impedance Z = 1/Y of the sample. The latter quantity reads, in terms of the reduced 
frequency y of eq. (2.7), 

^~ VCa;oy + 2iV2Ca;o^^"(y-y.)2 + l/(4Q2)- ^'-"'^ 

The real (dissipative) part of the impedance exhibits narrow resonance peaks, having 
Lorentzian line shapes, with a common absolute width 

Both the maximal value at resonance 

and the area under the resonance peak 

^"^L-. ^^^"^"^^ 2^ ^^-^"^ 

depend on the clusters and on the resonance under consideration only through the cross- 
section 7a. 

2.5. Duality symmetry 

Duality is one of the key concepts of the theory of planar graphs (see ref. [11] for 
an overview). It has far-reaching consequences in two-dimensional statistical physics (see 
ref. [12] for a review). The applications of duality to random resistor networks have 
been reviewed in ref. [1]. The dual of the network of Figure 1 has its electrodes along 
the horizontal sides (M and N are interchanged), and the dual clusters are obtained by 
drawing a dual bond across each bond of the original ones. Figure 2 shows some illustrative 
examples of clusters with their duals. Case (c) demonstrates that the dual of a connected 
cluster may be disconnected. 

The duality symmetry of the square lattice implies the identity 

Y{h) Y{h) = 1 (2.31) 

between the conductances of the original network and of its dual, with bond conductances 
related through ao = I/ctq, ?i = 1/cri, so that h = 1/h. By inserting the result (2.24) into 
the identity (2.31), we obtain the following predictions. The resonances of the dual set of 
clusters are located at the dual positions, 

ha = f/ha, i.e., Aa = 1 - Xa, Or UJa = ^\l^a, (2.32) 
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and any pair of dual resonances have identical cross-sections, 

la = la- (2.33) 

As an illustration we close up this general section by giving in Table 1 the positions 
Ao of the resonances, and the associated cross-sections 7a, for the three clusters shown in 
Figure 2. Cluster (a) is a generic example. Cluster (b) is self-dual, i.e., isometric to its 
dual, so that its resonances come in dual pairs. Cluster (c) has two peculiarities, namely 
it possesses a closed loop and its dual is not connected. 

3. APPLICATIONS 
3.1. One bond 

Consider first the simplest of all clusters, consisting of only one bond, joining the origin 
X = to x' = e^. The bond is parallel to the applied field for = 1, and perpendicular to 
it for /i = 2. The corresponding matrix M only involves two values of the Green's function. 
Go and Gi, given by eq. (A. 5). We thus get 




with eigenvalues 

Ai = 1/2, A2 = 0. (3.2) 
Only Ai leads to a resonance, lying at 

hi = -1, (3.3) 

i.e., Lo = u>o for the RL — C model. The associated cross-section reads 

71=45^,1' (3-4) 

with 5^ 1 being the Kronecker symbol. We thus have 71 = 4 for /i = 1 (the bond is parallel 
to the applied field), and 71 = for = 2 (the bond lies along an equipotential line). 

The latter property is quite general. Clusters consisting only of vertical bonds, perpen- 
dicular to the applied field, do not perturb it at all. Their resonances thus have vanishing 
cross-sections. 

3.2. Two bonds 

Consider now two bonds, in arbitrary relative position and orientation, namely one 
joining the points x = and x' = e^, and one joining the points y = (2/1,2/2) and 
y' = y -I- ej,, with = 1 or 2. 
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The corresponding matrix M is 



M 



/ -1/4 G^,y - G^y G^y - G^,y ^^ 

-1/4 1/4 Gx',y - Gx^y' Gx'.y' - Gx'.y 

G'x,y-G'x',y Gx'.y-Gx.y 1/4 "1/4 

VGx,y' — G^x',y' Gx',y' — Gx,y' —1/4 1/4 / 



(3.5) 



Its eigenvalues read 

Xi = l/2 + g, X2 = l/2-g, A3 = A4 = 0, (3.6) 

with 



9 — ^X,y — G-x^y' — Gy^/ y + Gx',y' 

= G'^,.(y) = G(y) - G(y + e,) - G{y - e^) + G{y + - e^). 
Only Ai and A2 lead to resonances, lying at 



(3.7) 



l-2g l + 2g 1 

= - lT2^' = -l^g = V,- ^'-'^ 

The associated cross-sections are 

2_iS^^^^ + S^ 2{5,,,-5.,,f 
^1- i_4^2 ' 72- . (3.9) 

In the limit where both bonds are infinitely far apart ((7 = 0), we obtain two degenerate 
resonances at /i = — 1. If the bonds are at a large but finite distance (|y| 3> 1), the estimate 
(A. 6) shows that the coupling between both resonances is approximately 

9 - f In |y| = -f - (5,,. - 2^^) « 1. (3.10) 

27Tdyf,dyt, 27r|y|^ V IyI / 

The resonances are thus symmetrically shifted by a small amount, 

hi,h2^ -l±4:g, (3.11) 
hence, in the case of the RL — C model, 

ui,LJ2 ^ u)o{lT2g). (3.12) 

3.3. Linear clusters 

Another case of interest is that of linear clusters. Consider a horizontal linear cluster 
consisting of ub = n consecutive bonds, joining the sites xei, with x = 0,---,n. The 
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entries of the corresponding matrix M only involve the Green's function G^, with the 
notation of the Appendix. They are symmetric for generic values of the site labels, 

M.x,y = 2Gx-y-Gx-y-i-Gx-y+i = F{x-y) (y = 1,- • -,71,- 1), (3.13) 

together with the non-symmetric boundary values 

^x,Q = Gx — Gx-li ^x,n = Gx-n — Gx-\-l-n- (3-14) 

The positions of the resonances and the corresponding cross-sections can be deter- 
mined analytically for the first few values of the size n: 

• n = 1. This is the one-bond case, investigated in section 3.1. 

• n — 2. This is a special case of the two-bond case investigated in section 3.2, correspond- 
ing to ji = V =1 and y = ei, hence 

£? = 2Gi-G2 = l/2-2/7r. (3.15) 

The general results (3.6, 8, 9) yield 

Ai = 1 - 2/7r = 0.36338, A2 = 2/7r = 0.63662, A3 = 0, 

/ii = -2/(7r - 2) = -1.7519, /i2 = 1 - 7r/2 = -0.57080, (3.16) 

71 = 7rV(7r - 2) = 8.6455, 72 = 0. 



n = 3. The matrix 



M = 



/ —Gi 2Gi — G2 2G2 — Gi — G3 G3 — G2 

G\ — 2G\ 2G\ — G2 G2 — Gx 

G2 — Gi 2G\ — G2 — 2G\ G\ 

\ G3 — G2 2G2 — Gi — G3 2Gi — G2 ~Gi 



(3.17) 



(3.18) 



has eigenvalues 

Ai = (7 - 16/7r -w)/4 = 0.28216, A2 = 8/7r - 2 = 0.54648, 
A3 = (7- 16/7r + t(;)/4 = 0.67136, A4 = 0, 

with w = ^384/7r2 - 224/7r + 33 = 0.77841, hence 

hi = -(16/7r + w- 3)/(7 - 16/7r - w) = -2.54411, 

h2 = -(3 - 8/7r)/(8/7r - 1) = -0.82990, (3.19) 

/i3 = -(16/7r - -u; - 3)/(7 - 16/7r + w) = -0.48951. 

We have also evaluated the associated cross-sections analytically, namely 

64(l + «;)(l-3/7r) 

71 = — ; ; ; — —r^ ; r = 14.633, 

' t(;(16/7r + t(;-7)(24/7r-t(;-7)(16/7r + t(;-3) 

72 = 0, (3.20) 

64(1 -w)(l- S/tt) 
~ w{-16/7T + W + 7)(24/7r + w- 7)(16/7r - w - 3) ~ ' 
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• n 3> 1. The resonance spectrum of very long linear clusters can also be investigated 
as follows. The bulk matrix entries (3.13) only involve an even function F{x — y) of the 
distance between nodes along the cluster. As a consequence, the matrix M is approximately 
a Toeplitz matrix, which is diagonalisable by means of the Fourier transformation. The 
eigenvalue A is given in terms of a wavevector (or momentum) q, by a dispersion relation 
of the form 

+00 

X(q)= J2 F{x)e-'i''. (3.21) 

x= — oo 

The expression (A. 2) of the Green's function leads to 

A(^) = r ^"'"'^ = = J'^^^^/^^' (3.22) 

'l + sm^{q/2) 



27r 2 — cos q — cosp V 3 — cos q 



For a large but finite linear cluster, the wavevector q is quantised. It assumes n 
discrete values in the range < g < tt, approximately given by 

nqpdaiT (a=l,---,n). (3.23) 

The positions of the resonances are therefore asymptotically distributed according to 
the smooth density 

p(\) = - -1 = (3.24) 

' TT dA 7r(l-A2)Vl-2A2 ^ ' 

Two limiting situations deserve our attention 

★ In the long- wavelength regime {q 0), which corresponds to low frequencies in the 
RL — C model, we have the linear dispersion law \{q) ~ q/2, hence the constant density 
of resonances 

p(A) = - (A ^ 0). (3.25) 

TT 

The scaling behaviour of the resonance spectrum in this regime will be investigated from 
a more general viewpoint in section 3.6. 

•k In the opposite limit (q — >^ tt), A approaches a non-trivial maximal value 

Al = l/v^ = 0.70711. (3.26) 

Near this upper band-edge of the dispersion law, the density of resonances (3.24) exhibits 
an inverse-square-root van-Hove divergency. 

The resonances of linear clusters thus extend over the range < A < A_l, i.e., —00 < 
h < Hl, with Hl = — (-\/2 — 1). The resonance spectra of the first 14 linear clusters are 
shown in Figure 3. 

We close up with a few words about the resonance cross-sections of linear clusters. 
First, the vanishing of 72 for the linear clusters with n = 2 and n = 3 [see respectively eqs. 
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(3.16) and (3.20)] is in fact quite general. Indeed, because of the left- right symmetry of 
linear clusters, the eigenvectors of the matrix M have a definite parity. We have Ra,n-k = 
(— l)"i?(j^fc, and a similar symmetry property for the L^, hence TZa = for a even. The 
cross-section 7^ therefore vanishes for every even a, and every n. Furthermore, in the 
long- wavelength regime (g <S 1, i.e., A <S 1 or a <S n), the eigenvectors La,k and Ra,k 
are given by superpositions of plane waves exp{±iqk) , and thus approximately uniformly 
extended over the cluster. Hence both amplitudes and TZa scale as the size n, so that 
we obtain the following scaling form of the cross-sections 

7a ^ r?Aa (a odd, n > 1), (3.27) 

where the Aa are numbers of order unity. This scaling law will be given a simple interpre- 
tation in section 4. 

3.4. Lattice animals 

A lattice animal (or an animal, for short) is by definition any connected cluster drawn 
on a given regular lattice. The problem of lattice animals has been an active field of 
statistical mechanics, one of the main motivations arising from cluster-expansion tech- 
niques for percolation and other problems in lattice statistics [13]. Just as percolation, the 
lattice-animal problem can be formulated either in terms of site occupation, or of bond 
occupation. The problem of bond animals is more suited for conduction properties. It con- 
sists in considering, with equal statistical weights, all clusters drawn on the lattice, made 
of a given number ub = n of bonds. The questions which have attracted most attention 
so far concern the asymptotic (n 00) behaviour of the total number of animals, and of 
geometrical characteristics, such as their mean radius of gyration. 

We have investigated the resonance spectra of animals drawn on the square lattice. 
To do so, we have generated all lattice animals, up to rimax = H, by adapting to the bond 
problem [10] an algorithm known for the site problem [14], and evaluated numerically the 
resonances of each animal. The total numbers of animals of size n are listed in Table 
2. This quantity grows asymptotically as 

it" 

iVn ~ — , (3.28) 

n 

where the universal 1/n power-law is an exact result in two dimensions [15], while the 
connectivity constant, here fx ~ 5.208, depends on the underlying lattice. 

Figure 4 shows histogram plots of the resonance spectra of all animals consisting of 
7, 9, and 11 bonds. Several peaks are clearly visible. Two of the most salient ones are 
shown on the plots, namely the upper band-edge Ay (3.26) of the linear clusters, and the 
upper band-edge of the T-shapcd clusters (i.e., the long linear clusters decorated by one 
bond branching perpendicularly at any place), that we have evaluated numerically to be 
at At = 0.83781. 
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3.5. Connections with the random binary model 



We now come to the question addressed in section 1, on the relationship between the 
singularities of the conductivity of the random binary network in the complex /i-plane, and 
the resonances of isolated clusters. We restrict again the analysis to the case of the square 
lattice. 

A systematic way of addressing this question is to consider the regime of very diluted 
impurities, 

c=l-p<l. (3.29) 

The conductivity has been shown [16-18] to admit a power-series expansion in c, of the 
form 

S(/i,c) = (7o(l + 6i(/i)c+62(/i)c2 + ---), (3.30) 

where the coefficients bk{h) only depend on the lattice under consideration and on the 
conductance ratio h. It is worth noticing that the series expansion (3.30) is different 
from the weak-disorder expansion of the conductivity as a power series in the successive 
cumulants of an arbitrary (smooth) distribution of the bond conductances, investigated in 
ref. [19]. 

The coefficients of the power-series expansion (3.30) are constrained by the duality 
symmetry of the square lattice. Indeed, by inserting this expansion into the thermodynamic 
version of the identity (2.31) [1], 

S(/i,c)E(l//i,c) = (j^, (3.31) 

independently of the concentration c of impurity bonds, we can derive an infinity of rela- 
tions for the functions bk{h), namely 

h{h) + bi{l/h) = 0, b2{h) + b2{l/h) = bi{hf, (3.32) 

and so on. 

• The first coefficient bi{h) only involves one-impurity effects. It is therefore correctly 
predicted by the EM A formula (1.6), which yields 

6,W=6r'^'(ft) = H^. (3.33) 

• The higher-order coefficients of the expansion (3.30) are far more difficult to calculate. 
The determination of b2{h), which contains two-impurity effects, has been performed on 
the cubic lattice [16], on the square lattice [17], and on the hypercubic lattice in any 
dimension D, by means of a two-impurity improved EMA-like scheme [18]. Coming back 
to the square lattice, the most appropriate expression of 62 (^) for the present purpose 
reads [18] 

b2{h)^h^{hr+b^{hf V — — 2 - 

2 (x,M)^o,i) 1 - Hh)G,,,{^)]' 1 - [b,{h)G,,,{^)]' 

(3.34) 
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This result is obviously more intricate than the EMA prediction 



Ah{h-lf _ h{hy{bi{h)+2) 
(/i + l)3 ~ 4 



(3.35) 



For instance, in the case of insulating impurity bonds {h = 0), eq. (3.34) yields 62(0) = 
-0.210749, while b^^^^\o) vanishes. 

First of all, the results (3.33, 34) obey the identities (3.32). This is obvious for 61, 
while more involved symmetry arguments are needed to check this property for b2, along 
the lines of ref. [19]. 

The result (3.33) shows that the coefficient bi{h) has a pole ai h = —1, corresponding 
to the resonance (3.3) of the one-bond problem. Similarly, eq. (3.34) shows that 62 (^) has 
an infinity of poles, situated at bi{h) = ±1/G^^j^(x), corresponding to all resonances of 
the two-bond configurations, determined in section 3.2. Indeed eq. (3.8) can be recast as 



This observation can be generalised to the following rule [10, 18]. For any k > 2, 
the coefficient bk{h) of the power-series expansion (3.30) has an infinity of poles along 
the negative real /i-axis, corresponding to all resonances of all sets of clusters consisting 
altogether of k bonds. This shows in particular that the conductivity of the random binary 
model is singular along the whole negative real axis of the complex /i-variable. Indeed our 
investigations of the resonance spectra of linear clusters and of lattice animals demonstrate 
a dense accumulation on the negative real axis. 

3.6. Fractal clusters 

This last section is devoted to the resonances of large clusters drawn on the square lat- 
tice, including fractal ones. It turns out that the spectra of large fractal clusters generically 
exhibit a scaling power-law behaviour in the A <^ 1 regime, i.e., for a -C ur. We shall ex- 
plain this phenomenon in a heuristic way, making use of our intuition of the low-frequency 
response of RL — C clusters, where A ^ {uj/ujof'. 

Consider a large but finite patch of a fractal cluster, of diameter £, embedded in the 
square lattice. The numbers of its sites, bonds, and resonances scale as 



where Di? < 2 is the Hausdorff (or fractal) dimension of the cluster. 

The investigation of the resonances of linear clusters, performed in section 3.3, has 
shown that the resonances for A <S 1 are characterised by coherent long-wavelength ex- 
citations along the cluster. We propose to generalise this picture to clusters of arbitrary 
shape, and to assert that the lowest resonance is generically given by 



b^{h) = ±l/g. 



(3.36) 



A 



(<^min/^o) 



1 



(3.37) 



'mm 



|2' 
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where L{£) and C{£) are the effective inductance and the effective capacitance of the 
cluster, considered as a whole: 

• The effective inductance L{£) of a cluster obeys the same scaling behaviour as its end- 
to-end resistance, 

L{£) ~ R{£) ~ f/", (3.38) 

where t/u is the usual notation coming from finite-size scaling theory for the percolation 
problem (see e.g. ref. [1]). 

• To the contrary, the capacitance C{£) of the cluster is assumed to be insensitive to its 
internal structure. In two dimensions it is therefore independent of £, so that we obtain 

Amin ~ r'/". (3.39) 

This scaling law can be extended to all resonances such that A <^ 1, i.e., a <^ nji, by 
expressing that in this regime A^ only depends on the label a through the dimensionless 
combination a/riR -C 1. We thus get 

Aa-(a/nfi)^ (a<nK), with ( = . (3.40) 

u Up 

An equivalent statement consists in writing the scaling law 

p(A) - A-^+^/^ (A ^ 0) (3.41) 

for the density of resonances of the infinite fractal structure. 

For linear clusters, we have t/u — Dp — 1, so that C = 1, in agreement with the 
analytical results (3.23, 24). 

We have also investigated three examples of deterministic self-similar fractal clusters, 
which we call the Worm, the Cross and the Gasket. Their iterative construction rules are 
shown in Figure 5, while the resulting clusters are shown in Figure 6. At each step of the 
construction, the number of bonds is increased by a constant factor a, while the diameter 
of the cluster (in units of the bond length) is multiplied by a scaling factor b. After k steps 
of iteration, the fractal cluster of the A;-th generation thus consists of ub ~ bonds, while 
its diameter scales as £ ~ 6^^. As a consequence, the fractal dimension of these clusters is 

= lna/ln6. (3.42) 

The exponent t/v of the end-to-end resistance is trivial in the first two examples, which 
possess no closed loops. Its value for the Gasket is known for long, and rederived e.g. 
in the review [1]. The values of the various scaling factors and exponents for the three 
fractal clusters are listed in Table 3. We have evaluated numerically the resonance spectra 
of successive generations of these three fractal clusters, by diagonalising the associated M 
matrices, up to a maximal scale A;max corresponding to a few thousand bonds. Figure 7 
shows logarithmic plots of the resonance positions against a/uR <^ 1, for the largest two 
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scales of each type of cluster. Power laws clearly show up, with exponents C in agreement 
with the prediction (3.40), listed in Table 3. 

4. DISCUSSION 

We have investigated the resonances which show up in the conductivity of a connected 
cluster, and of a finite set of clusters, drawn on the square lattice and embedded in a large 
rectangular sample. Our central result (2.24) shows that the conductivity is characterised 
by a finite number ur of resonances, located at well-defined negative real values ha of the 
ratio h of the conductivities of both phases. This formalism is chiefly aimed at describing 
the resonant dielectric response of clusters within the RL — C model. 

The strength of each resonance is measured by a cross-section 7a, given by eq. (2.25). 
The cross-section is naturally interpreted as the area of the sample over which the cluster 
significantly perturbs the applied field £ at resonance. A slightly different and comple- 
mentary viewpoint is as follows. Each factor Ca or Ra of eq. (2.25) has the dimension of 
a length. It can be interpreted as a measure of the strength of the dipole induced on the 
cluster by the applied field, so that the cross-section 7^ is proportional to the square of the 
induced dipole. The result (3.20) shows however that there is no simple algebraic relation 
between Ca and TZa in general, and that 7^ is not mathematically a perfect square. 

This interpretation of the cross-section sheds some light on several aspects of this work. 
First, the coupling constant g (3.10) between the resonances of two distant bonds coincides 
with the interaction energy of two dipoles D and D' in two-dimensional electrostatics, 
namely 

. = -,0.V),0-.V')^ = 4(^-.(^). (4.) 

The n^-law (3.27) for the cross-sections of long linear clusters is also simply understood 
as follows. For A -C 1, i.e., at low frequency for the RL — C model, the clusters respond 
coherently, so that the induced dipole is proportional to the size n. Our interpretation 
of the cross-section also allows a qualitative discussion of the finite-size corrections to 
our predictions. The most important of these effects consists in a shift of the spectrum 
of resonances. For a large but finite sample of size M x N, the elements of the matrix 
M differ from those corresponding to the infinite square lattice by terms of order 1/M^ 
and I/N'^i proportional to the interaction energy (4.1) between the dipole induced on the 
cluster and its mirror images with respect to the boundaries of the sample. The resonances 
are therefore shifted by a small amount of order 1/M^ and with respect to their 

theoretical positions ha- 

The present analysis can be extended to the problem of a coloured cluster, consisting of 
different kinds of metallic bonds. Assume that each bond of the cluster has a conductance 
c"x,y = c"o/ix,y, where the /ix,y are given arbitrary complex numbers. A generalisation of 
eq. (2.10) shows that the resonance condition now reads 

det(l-N) = 0, (4.2) 
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where N is an ns x ns matrix, defined by 



N. 



^ ^y.z) (^x,y - Gx,z) • 



(4.3) 



z€C(y) 



The condition (4.2) yields a polynomial relation between the /ix,y A simple illustration is 
given by the linear cluster consisting of two neighbouring bonds with conductance ratios 
hi and /i2. We thus obtain the relation 



If all bonds of the cluster are identical, i.e., /ix,y = h, we have N = M/A. The eigenvalue 
criterion is thus recovered. In our example, eq. (4.4) for hi = /12 yields back the result 



We have also investigated the small-A region of the resonance spectra of large clusters, 
including fractal ones, embedded in the square lattice. The scaling laws (3.40), (3.41), es- 
tablished in a heuristic way, have been checked analytically by means of the dispersion law 
pertaining to linear clusters, and numerically on several examples of self-similar fractals. 
We would like to emphasise that these scaling laws cannot be generalised in a straightfor- 
ward way to clusters embedded in a higher-dimensional lattice. This can be demonstrated 
explicitly in the case of linear clusters. On the hypercubic lattice in dimension D, the 
dispersion relation for a small wavevector {q <^ 1) reads 



This result implies the existence of an upper critical dimension Dc — 3, at least in the case 
of linear clusters. 

It should also be noticed that the present results only concern the dielectric resonances 
of finite clusters embedded in an infinite matrix, which consists of a regular lattice. The 
resonance spectra of deterministic models for fractal structures, such as the incipient per- 
colation cluster, investigated in refs. [20, 21], have very different characteristics. In the 
latter case both phases have comparable sizes and geometries, so that the notions of clus- 
ter and matrix are no more pertinent. The spectra of such model systems are generically 
themselves fractal. They are supported by a Cantor set of the negative real /i-axis, instead 
of the whole axis in the present situation. 

To close up, we come back to our initial motivation, namely a better understanding 
of the analytic structure of the conductivity of the random binary model, both in the 
conductance ratio h and in the concentration variable p or c. First of all, and from a 
qualitative viewpoint, the present analysis confirms that the cut of the conductivity in 
the complex /i-plane can be viewed as the accumulation of the resonances of clusters and 
sets of clusters with any size and shape, and that it extends over the whole negative real 
axis. The more quantitative analysis of section 3.5 shows that every coefficient bk{h) of 



(1 + hi){l + h2) = (1 - 4/7r)2(l - hi){l - h2). 



(4.4) 



(3.16). 




(4.5) 
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the expansion (3.30) of the conductivity has itself a countable infinity of poles, associated 
with the resonances of all embcddings, connected or not, of k bonds in the square lattice. 
The present work leaves several open questions, such as e.g. a more accurate form for the 
Lifshitz tail (1.9). It however shows in a very suggestive way how intricate the exact con- 
ductivity of the binary model can be on a finite-dimensional lattice, especially as compared 
to the simple and very commonly used EMA formula. 
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APPENDIX 

In this Appendix we summarise the main properties of the Green's function of the 
finite-difference Laplace operator A on the square lattice, which are useful in the body of 
this paper. The reader is referred to ref. [22] for a more comprehensive exposition. 

The Green's function Gx,y = G(x — y) is by definition a solution of 

-AG(x) = S^,o, (A.l) 

with 5x,o being the Kronecker symbol. The difference equation (A.l) has a unique solution 
with all required symmetries, up to an additive constant, which we fix by setting G{0) = 0. 
We obtain by Fourier transformation 

f d^p e^P-^''"'') - 1 
Gx,y = ^(x - y) = G{xi - yi, X2 - 2/2) = , (A.2) 

where the double integral runs over the first Brillouin zone B (— tt < pi,P2 < tt), and 
where 

K(p) = 2(2 - cospi - C0SP2) (A.3) 

is the Fourier transform of (—A). 

The values of (j(x) along the diagonals, namely for m = ±n, can be evaluated explic- 
itly from the representation (A.2) by means of elementary integrals. One thus obtains 

G(±n, ±n) = (^1 + i + ■ ■ ■ + ^-^^ (n > 1). (A.4) 

The values of G(x) all over the square lattice can then be determined recursively from eq. 
(A.l), starting from the knowledge of the values (A.4) [22]. 

Introducing the short-hand notation Gm — G{±m, 0) = G{0, ±m) for the values of 
the Green's function along the co-ordinate axes, we have 

Go = 0, Gi = -1/4, G2 = 2/7r - 1 = -0.36338, G3 = 12/7r - 17/4 = -0.43028, 

(A.5) 

and so on. 

Finally, the Green's function admits the long-distance expansion 



G(x) - — 

^ ' 27r 



ln|x| + (-ln2 + 7. + — + 



(|x| > 1). (A.6) 



The leading isotropic logarithmic term is nothing but the Green's function of the differential 
Laplace operator in the plane, which represents e.g. the potential of a point charge in two- 
dimensional electrostatics. The finite part involves Euler's constant 7^ = 0.57721. The 
first anisotropic correction, due to lattice effects, is of relative order l/|xp. 
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CAPTIONS OF FIGURES AND TABLES 



Figure 1: Schema of the sample used in this work. 

Figure 2: Three examples of clusters (full lines), together with their duals (dotted lines). 

Figure 3: Plot of the resonance positions of the first 14 linear clusters. The dashed line 
shows the upper band-edge Al of eq. (3.26). 

Figure 4: Histogram plots of the resonance spectra of all lattice animals with (a) 7 bonds, 
(b) 9 bonds, and (c) 11 bonds. The letters show the upper band-edges Xl and At of linear 
and T-shaped clusters, respectively. 

Figure 5: Construction rules of the fractal clusters considered in section 3.6. 

Figure 6: Largest generation whose resonance spectrum has been evaluated numerically, 
for (a) the Worm, (b) the Cross, and (c) the Gasket. 

Figure 7: Logarithmic plot of the small-A range of the resonance spectra of the largest 
two generations of (a) the Worm, (b) the Cross, and (c) the Gasket. The slopes of the full 
lines are the exact theoretical values of listed in Table 3. 

Table 1: Positions of the resonances Ao, and associated cross-sections 7^, of the clusters 
shown in Figure 2. 

Table 2: Total numbers of lattice animals of n bonds, up to ?imax = 11. 

Table 3: Scaling factors and exponents of the fractal clusters considered in section 3.6. 
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Table 1 



cluster 


resonances 


cross-sections 


(a) 


Ai = 0.23605 
A2 = 0.31285 
A3 = 0.59706 
A4 = 0.85404 


71 = 5.18968 

72 = 2.17581 

73 = 2.42350 

74 = 0.10714 


(b) 


Ai = 0.13067 
A2 = 0.30756 
A3 = 0.50000 
A4 = 0.69244 
A5 = 0.86933 


71 = 1.08684 

72 = 0.48889 

73 = 10.1794 

74 = 0.48889 

75 = 1.08684 


(c) 


Ai = 0.24692 
A2 = 0.58281 

A3 = 0.64009 
A4 = 0.68901 
A5 = 0.84116 


71 = 12.9772 

72 = 1.11399 

73 = 5.17747 

74 = 0.49606 

75 = 0.12694 
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Table 2 



71 




1 


2 


2 


6 


3 


22 


4 


88 


5 


372 


6 


1628 


7 


7312 


8 


33466 


9 


155446 


10 


730 534 


11 


3466170 



Table 3 



Fractal 


a 


b 


Dp 


t/v 


c 


Worm 


3 




21n3/ln5 
= 1.36521 


21n3/ ln5 
= 1.36521 


1 


Cross 


5 


3 


In5/ln3 
= 1.46497 


1 


In3/ln5 
= 0.68261 


Gasket 


3 


2 


In 3/ In 2 
= 1.58496 


ln(5/3)/ln2 
= 0.73697 


ln(5/3)/ln3 
= 0.46497 
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